ON A COUPLED PDE MODEL FOR IMAGE 
RESTORATION 



V.B. SURYA PRASATH AND DMITRY VOROTNIKOV 

Abstract. In this paper, we consider a new coupled PDE model 
for image restoration. Both the image and the edge variables are 
incorporated by coupling them into two different PDEs. It is shown 
that the initial-boundary value problem has global in time dissi- 
pative solutions (in a sense going back to P.-L. Lions), and several 
properties of these solutions are established. This is a rough draft, 
and the final version of the paper will contain a modelling part and 
numerical experiments. 



1. Introduction 
We consider the following problem 

(1) ^^^ = div{g{v{t,x))Vu{t,x)), 

(2) - \ix)Avit,x) = (1 - A(x))(|Vu(t,x)| - vit,x)), 



(3) u 



0,v 



an ' 



= 0, 

an 



(4) u\t=o = Uo, v\t=o = Vq. 

It is shown that the problem possesses global in time dissipative 
solutions] uniqueness, regularity and some other properties of these 
solutions are studied. The concept of dissipative solution was suggested 
in [7] for the Euler equations of ideal fluid flow, which are still not 
proven to have global weak solvability. Later, existence of dissipative 
solutions was established for Boltzmann's equation [6l H], the ideal 
MHD equations [12], Navier-Stokes-Maxwell equations [2], Euler-a and 
Maxwell-a models [10] and viscoelastic diffusion equations [TT] . 



Key words and phrases, global solvability, dissipative solution, imaging. 
Research partially supported by CMUC and FCT (Portugal), through European 
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The features of our problem ([Il)-(jll) which oppose strong and classical 
weak well-posedness are the presence of a nonlinear function (modu- 
lus) of the gradient of u in the right-hand side of ([2]) and the Perona- 
Malik-like form of g. The inequality (12. ip in the definition of dissipa- 
tive solutions turns out to contain the absolute value function as well. 
Therefore, unlike in the previous works on dissipative solutions, it is 
impossible to pass to the limit in this inequality via weak and weak-* 
compactness argument. Nevertheless, we manage to do it via strong 
compactness, although it is not sufficiently strong to obtain classical 
(i.e. not dissipative) weak solutions. 



2. Well-posedness of the problem 

The objective of this section is to prove Theorem 12.11 concerning ex- 
istence, uniqueness, regularity and some other properties of dissipative 
solutions to the problem ([I])-(j4]). We consider the simplest Dirichlet 
boundary condition ([3]), but other boundary conditions can also be 
handled. 

In the section, Q is considered to be a domain (i.e. an open set in 
M^) possessing the cone property. We recall pQ that this means that 
each point x & Q is a vertex of a finite cone contained in Q, and all 
these cones Cx are congruent. A finite cone is a set of the form 

C, = Bin{x + ^{y -x)\ye B^, ^ > 0} 

where Bi and B2 are open balls in R^, Bi is centered at x, and B2 does 
not contain x. Obviously, rectangular domains have this property. 

The symbol C will stand for a generic positive constant that can 
take different values in different lines. We sometimes write C{. . .) to 
specify that the constant depends on a certain parameter or value. 

We assume that (7 : R — )■ R, -i= and A : i7 — ^ R are Lipschitz 
functions having positive values, g is bounded, A < 1, 

(5) Ao = inf \{x) > 0. 

The assumptions on g hold, for instance, if 



(6) g{s) 



a 



6 + c|s|'^' 



where a, b, c, d are positive numbers, and 1 < d < 2. 
Note that 



1 



1 



V9{s) V9{s) ^/g{0) V^(0) 
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We use the standard notations Lp(fi), W^{n), = W^{n) 

for the Lebesgue and Sobolev spaces. We will often keep the function 
space symbol and omit Q. 

The Euclidean norm in finite-dimensional spaces is denoted by | ■ |. 
The symbol || ■ || will stand for the Euclidean norm in L2{fl). The 
corresponding scalar products is denoted by a dot ■ and parentheses 

(-,■)• 

Let Hq(Q) be the closure of the set of smooth, compactly supported 
in Q, functions in H^{Q). By virtue of Friedrichs' inequality, the Eu- 
clidean norm || ■ ||i corresponding to the scalar product 

(m, f )i = (Vm, Vf ) 

is a norm in Hq. 

The set V2 = H^^fl) fl H'^{Q) is a Hilbert space with the scalar 
product 

\a\=2 

Denote the corresponding Euclidean norm by || ■ II2. 
Let Vr, I < r < 2, be the closure of V2 in W^. 

We recall the following abstract observation [9l [13] . Assume that we 
have two Hilbert spaces, X C Y, with continuous embedding operator 
i : X ^ Y, and i{X) is dense in Y. The adjoint operator i* : Y* ^ X* 
is continuous and, since i{X) is dense in Y, one-to-one. Since i is one- 
to-one, i*{Y*) is dense in X*, and one may identify Y* with a dense 
subspace of X*. Due to the Riesz representation theorem, one may 
also identify Y with Y*. We arrive at the chain of inclusions: 

(8) X CY = Y* C X*. 

Both embeddings here are dense and continuous. Observe that in this 
situation, for / G F, m G X, their scalar product in Y coincides with 
the value of the functional / from X* on the element m G X: 

(9) (/,M)y = (/,n). 

Such triples (X, F, X*) are called Lions triples. We use the Lions triples 
iV^^L^.V*) and {HlL^.R-^). 

The symbols C{J] E), C^^J'; E), L2{J; E) etc. denote the spaces of 
continuous, weakly continuous, quadratically integrable etc. functions 
on an interval JT" C M with values in a Banach space E. We recall that 
a function u : J ^ E is weakly continuous if for any linear continuous 
functional g on E the function g{u{-)) : J' ^ M. is continuous. 

We require the following spaces 

= W,{Q,T) = {re L2{0,T;V2), r' G 1.2(0, r;^*)}. 
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Ikllwi = ||T||L2(0,T;y2) + l|T'||L2(0,r;V2*); 

W2 = W2{^,T) = {r G L2{Q,T-Hl), r' G L^iQ.T-R-')}, 

Ikl|w2 = \\A\l2{0,T;HI) + lk'||L2(0,T;H-i)- 

Let us introduce the operator 

where ip is an arbitrary element of V2. 

Denote by TZ the following class of pairs of functions: 

n = u,ioc{o, 00; V2) n Loo(o, 00; w^) n w^^M oo- L2) 

xL2,ioc(0, 00; V2) n Loo(0, OO] Loo) n W2^,;oc(0, oo; L2). 
Observe that the following expressions, where 5 is a positive number, 
are well-defined for (w, r) G 7^, and their values are in L2,ioc(0, oo; L2): 

(10) i?i(«;,r,5) = + 5div((7(r)Vw;), 

(11) E2{w,T,5) = -^ + AAr + 5(l-A)(|Vt/7|-r) + (l-5)(Vr-VA), 

(12) Ei{w,t) = E^{w,T,l), 

(13) E2(w,r) = E2(^i;,r,l). 
Let us recall the Sobolev inequality 

(14) \\u\\l^<C(SI)\M2. ueV2. 
and the Ladyzhenskaya inequahty [13] 

(15) \\n^\\ < \/2\\u\\ ||Vm||, u G i^o^. 

The following Gronwall-like lemma will be useful. 

Lemma 2.1. Lemma 3.1]; Let f,x,L,M : [0,T] ^ R be scalar 
functions, x,L,M G Li(0,T), anc? / G Wi{0,T) (i.e. f is absolutely 
continuous). If 

(16) x(t)>0,L(t)>0 
and 

(17) f{t) + x{t)<mf{t) + M{t) 
for almost all t G (0, T), then 

t 

f{t)+ [ xis)ds< 
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/(0) + J^wlj M{s)ds 



(18) exp I j L{s)ds 

.0 

for allte [0,T]. 
We can now give 

Definition 2.1. Let Uo,Vq G L2{^1). A pair of functions {u,v) from 
the class 

(19) u,v eC^{[0,ooy,L2), 

is called a dissipative solution to problem dH) - (HI) if , for all test func- 
tions {(, 9) and all non-negative moments of time t, one has 

7H*)ii^ [\Ht)-m? + \Ht)-e{t)r] 



<7'*+""""'{lko-C(0)f + ||t;o 



(20) + 2r' {E,{c.e){s)Ms)-as)) + [E2{c.e){s),v{s)-e{s)) 



where 7= 7(fi, g, A, C, 6') > 1 is a certain function of Vt, g, A, ( and 9. 

Theorem 2.1. a) Given uq,vq G L2, there is a dissipative solution to 
problem - 

b) This solution {u, v) belongs to L4/3 ;oc(0, oo; \^_e+4/3) XL2 ioc(0, 00; Hq), 

0<6<|. 

c) If, for some Uq^Vq G L2, there exist T > and a strong solution 
{ut,vt) to problem ([I]) - dl]), which is a restriction of a function from 
IZ to (0,T). Then the restriction of any dissipative solution (with the 
same initial data) to (0,T) coincides with {ut,Vt)- 

d) Every strong solution {u,v) & is a (unique) dissipative solution. 

e) The dissipative solutions satisfy the initial condition (jlj). 

To prove Theorem 12.11 we consider the following auxiliary problem: 



(21) 



— + eAu = S div{g{v)'Vu) 



dv 



(22) —-XAv = S{1 - A)(|Vm| -v) + {1~ 5){Vv ■ VA), 

L/t' 



(23) 



u 



on 



0,v 
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(24) u\t=o = Suo, v\t=o = 6vo. 

Here, e > and < 5 < 1 are parameters. 

The weak formulation of (12T|) - (l24l) is as follows. 

Definition 2.2. A pair of functions {u,v) from the class 

(25) ueWi, V eW2 

is a weak solution to problem fl2T]) - flM|) if the equalities 

(26) ^(m, (^) + £(m, ifi)2 + 5igiv)Vu, Vip) = 0, 
and 

(27) +5(V^;, 0VA) - 5 ((1 - A)(|Vn| - i;), 0) = 

are satisfied for all G V2, G //g almost everywhere in (0,T), and 
(El and (IHD hold. 



Lemma 2.2. Let {u,v) be a weak solution to problem ( !2T]) - ([24 

Then, for all test functions {(, 9) and < t < T , one has 

^Mtw^\\u{t)-cit)f + \Ht)-mf 

+2e J \\uis)-Cis)\\lds + Xo J \\v{s)-e{s)\\lds} 


< 72*+^ll"ol!^(||5«o - C(0)ir + l|5t;o - e(0)f 



t 



+ / 27 







(28) +(E2(C, 0, 5){s),v{s) - e{s)) - £(C(s), u{s) - C(s))2| ds] 
where 7= 7(fi, (7, A, C, ^^) > 1 zs a certain function of Vt, g, A, C anc? 0. 

Proof. Let us first derive the straightforward energy estimate. For al- 
most all t e (0,T), let ^ = u{t) in TherQ 

(29) ~{u, u) + 6igiv)Vu, Vu) + eiu, u)2 = 0. 



^See e.g. [131 p. 153] on how ^ appears in (f^ . 
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Integration in time gives 



t 



(30) -\\u{t)r+ {6g{v{s))Vu{s),Vu{s))ds<-\\uor. 







Observe now that 



(31) ^(C, ^) + 6{g{e)VC, V^) + {E,{C e, 6),ip)+ e{C, y.)^ = e(C, ¥^)2, 
and 

(32) -5 ((1 - A)(| vci -e)A) + mc. 0, 6),<P) = 0. 

for € V2, € Hq. Denote w = u — Q and q = v — 6. For almost all 
t G (0, T), put ip = w{t) and = ^"(t). Add the difference between fl26|) 
and fl3T]) with the difference between (12 ■2p and ([2]), arriving at 

^ ^^(^,^) + '^(^(^)Vw, Vw) 

+e(u;, u;)2 + (AV?, V^) + S {{1 - A)?, ^) 
= -6{[g{v) - g{e)]VC, Vw) + 5 ((1 - X){\Vu\ - | VC|), ^) - 5(V?, ^VA) 

(33) +(Ei(C, ^, 5), + (i?2(C, 5), ^) - £(C, ti^)2. 

Let us estimate the first three terms in the right-hand side. 
- 6{[g{v) - g{e)]VC, Vw) + 6 {{1 - X){\Vu\ - | VC|), ^) 
<C{C,g)6{\v-ei\Vw\) 

kl 



<CiC9){^=,VS9F)\'^w\ 



< C{C,e,g) fkl, VM^)\Vw\) + C{C,g) U\ ^/6^){\VC\ + \Vu\) 



(34) < II v/5^V«;f + CiC, e, g) |kf + C(C, ^5^1 Vw| 
and 

(35) -5(V^,,VA) < C(A)(^,Vc:) < ^||,||2 + c(A)|k||2 
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Now, ([2]) implies 

Id, . Id 3Ao II ||2 

<C{(,9,\,9) [e,l + VS^\Vu 

(36) +(Ei(C, 9, 6),w) + {E^iC, 9, 5), ^) - ^(C, w^. 

Denote <l>{t) = 1 + ^/6g{v{t))\Vu{t)\ . Due to 

f^iw,w) + |(,,,) + 2e{w,wh + ^llVc^f 
<C{C,9,X,g)m\\\V4 
+2(Ei(C, 9, S),w) + 2(E2(C, 0, 5),^)- 2£(C, ^)2. 



(37) 
Thus, 



2||^l|2 



(38) 



<C(C,^,A,^?)$1^ 

+2(Ei(C, ^, 5), w) + 2(E2(C, ^, 5), - 2£(C, ^)2. 

We now require two estimates for $, 

t 

<^\s)ds = j + ^5g{v{sy)\Vu{s)\fdxds 

Q 

<2y y dxds + 2 j j 5g{v{s))\Vu{s)\^ dxds 
on on 



(39) < 2t|(]| +(5||mo|| 
by virtue of (l30l) . and 

(40) j^\s)ds> j j dxds = t\n\. 
on 

With the help of Lemma 12.11 we derive from ([2])- pUj) that 

t i 

IkWir+lkWf + 2£ / |k(s)||^ds + Ao / \\V^is)fds 
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<exp (c{C,e,X,g) J <l>\s)ds\ |||^(0)f + ||,(0)||2+ 









exp \^C{(,9,X,g) j <l>\Od^ j [2{E^{(,9,5){s),w{s)) 

+2{E2{C, 9, S){s), ^{s)) - 2sias),w{s))2] ds] 
< exp {C{C,9,X,g){2t\n\ + 5||nof - )) {|h(0)f + |k(0)f + 



exp (-C(C, 9, A, g)s\Q\) \2{E^{C, 9, 6){s), w{s)) 



+2{E2{C, 9, 6){s),<;{s)) - 2e{C{s),w{s))2\ ds ^ 
< exp (C(C, 9, X,gm\ + l)(2t + ^^of - IkWf )) |lk(0)|P+|k(0)|p+ 



exp (-C(C, 9, A, g)si\n\ + 1)) |2(Ei(C, ^, ^(.)) 



(41) +2(E2(C, ^(s)) - 2Eias),wis))2\ rfsj, 

since s < 2t. Now ([2]) yields (ES]) with 

7 = exp{C(C,^,A,(7)(|fi| + l)}. 

□ 

Lemma 2.3. Let {u, v) be a weak solution to problem fl?I]) - fl2^ . T/ie 
following estimates are valid: 

(42) hl|Loo(0,T;L2) + II^IUoo(0,T;L2) + || ^ || L^CO.T;//!) < C, 

c 

(43) ||^^||L2(o,r;y2) < -7=^ 



(44) ||VM||i2(o,T;Li) + ||Vm||li(0,T;L,.) + \\^u\\ L^^^^o^T;L_,+^/3) < C, 



1 

l<r<2,0<e<-, 
o 



(45) lk'IU2(0,T;V2-) + \W\\l2{0,T;H-2) < (1 + V^)^: 

(46) ||t;'|U,(o,r;H-i) < (1 + l/v^)C^- 
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The constants C = C(T,\\uo\\,\\vo\\, X, g,^) are independent of e and 
S. 

Proof. The estimates fH2l) and fj43|) are direct consequences of (12.21) 
with C = ^ = 0. 

Then, using ([7]) and ( l30i) . we have 

||VM||i2(o,T;Li) < 
II V^9iv)^u\\L2{0,T;L2) II V IU^(0,T;L2) 

(47) <C\\1 + \v\\\lml.)<C, 

and, since ifg C Lp for any p < oo by Sobolev embedding, 

\\^u\\L,{0,T;Lr) < 
II V^5^Vm||l2(o,T;L2)||1 + |t'IIU2(0,T;L2,,/(2_,)) 

(48) <C(l + ||i;|L,(o,T;/^o^))<a 

By the time-space Holder inequahty [lU Lemma 2.2.1(b)], 

l|V«|U,/3(0,T;L_.,,/3)<ll|V«r/1 L4{0,T; L.)|||V^r/l 

(49) < J||Vn|U,(0,T;L,)l|V^i||L,(0,T;L4-3.) < C. 

V "3^3? 

It remains to estimate the time derivatives, expressing them from 
m\) and ([22D- Utihzing we get 

\\{u',^)\\L2io,T) < 5||(c/(t;)VM, V(/?)|U2(o,r) +e||(M,</')2||L2(o,r) 

(50) <C(l + v^)||^||2, 
and 

||(t^',</')||L2(o,T) < ||(AVt;,V0)|U,(o,T) + 5||(Vt;,0VA)|U,(o,T) 
((1 - A)t;,0) |U,(o,T) + 5|| ((1 - A)|VnU) |U,(o,t) 

< ll^llL2(O,r;/fi)ll0l|l + C(A)||t^||L2(O,T;i/i)ll0ll 

(51) +||V«|U,(o,T;L.)II0I|l^<C||0||2. 

In order to get ( H6i) . it suffices to observe that 

5||((l-A)|VwU)|U,(o,T) 

(52) < ||V«1U,(O,T;L2)||0|| < C\\uU.,iO,T;Vm\l < ^Ml- 

□ 
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Lemma 2.4. Given T > and uq, vq G L2, there exists a weak solution 
to problem ([2ID - dMD with 6 = 1. 

Proof. Let us rewrite the weak statement of (!2T]) - (12^ in the suitable 
operator form 

(53) A{u,v) = 5Q{u,v). 

The operators i, Q : W^i x W2 ^ ^2(0, T; V^) x ^2(0, T; H'^) x L2 x L2 
are determined by the formulas 

{A{u,v),{^,(f))) 
= (^("' + ^(^' '^)2' ^(^5 0) + (AVi;, V0), M|i=o, w|t=o) , 



{Q{u,v), (¥?,0)) 

= ( - {9iv)Vu, V^), -(Vt;, 0VA) + ((1 - A)(| Vu| -v),(j)), uq, vq^ . 

Here ip E V2 and (p G -ffg ^^^t functions. 

The operator Q is continuous and compact. Here we only explain 
this claim for its first component, and for the others the proof is 
more straightforward. We observe first that the embedding Wi C 
Lp{0,T;Wp) is compact for some p > 2. This can be shown using 
[H Corollary 8]. The embedding W2 C L2{0,T; L2) is compact by [H 
Corollary 4]. Let {um,Vm) {uo,vo) be a weakly converging sequence 
in Wi X W2- Then («„,!>„) is strongly converging in Lp{0,T;Wp) x 
L2{0,T; L2). By Krasnoselskii's theorem [51 Theorem 2.1], g{vm) — ^ 
g{vo) in Lg(0,T;Lg) for any q < +00. Thus, g{vm)Vum g{vo)VuQ 
in L2(0,T; L2), and the claim follows. 

The linear operator A is continuous by [13l Corollary 2.2.3] and 
invertible by [131 Lemma 3.1.3]. Thus, fl5^ can be rewritten as 

(54) {u,v) = 5A-^Q{u,v) 

in the space Wi x W2- 

Lemma 12.31 yields the a priori estimate 

(55) \\u\\wi + \\v\\w2 < 

where C may depend on e but does not depend on 5. By Schaeffer's 
theorem ^ p. 539], there exists a fixed point of the map A^^Q, which 
is the required solution. □ 



We will also need the following simple fact. 
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Proposition 2.1. Let G be a measurable set in a finite- dimensional 
space, X : M — 7- M 6e a continuous function, and let t/m '■ G ^ M. be 
a sequence of functions. Assume that {ym} is uniformly bounded in 
L^{G), and ym yo m Lg{G), q>l. Then xiVm) xivo) m Lp{G) 
for any p < oo. 

Proof. Due to the uniform boundedness of {ym}, without loss of gen- 
erality we may assume that x is also bounded, and then it suffices to 
apply [5i Theorem 2.1]. □ 

Based on the obtained lemmas, we can proceed with the sketch of 
the proof of Theorem 12.11 We refer to [10] for the details of the 
technique, and mainly focus on the new issues. To prove a) and b), 
one passes to the limit in (12. 2 p with 5 = 1 as e = on every 

interval (0, T), T > 0. However, unlike in [71 [121 IHl ID], in view of the 
presence of the absolute value in the right-hand member of f l2.2p . it is 
not possible to do it via weak and weak-* compactness. 

Let {um, Vm) be the weak solution to problem fl2T]) - with e = 
Em- Lemma 12. Sj [H Corollary 4] and the compact Sobolev embed- 
ding W^^^^^^ C L2 imply that without loss of generahty Um u in 
174/3(0, T; L2), Vm ^ V in L2{0,T; L2). Then, by fH2]) and Proposition 

in L2(0,T). Furthermore, by the same proposition, ||um(t) — C(^)ll^ ~^ 
\\u{t) - C{t)f, \\vm{t) - e{t)f -> \\v{t) - e{t)f in L2(0,T). Therefore 

^\\-r.m' {\\Umit) - + \\Vrnit) - 9{t)r} 

(56) ^ ^Mtr {\\u{t) - cm' + \Ht) - eit)f} 

in Li(0,T). Note that 

(57) e e L4(0, T; H') C 1^(0, T; L2) H ^2(0, T; H''). 

This yields El (C,^) G L^{0,T]L2). Remember that E2(C, ^) G L2(0,T;L2). 
Thus, we can pass to the limit in the right-hand side of (12.20 as well; 
the last summand (the one with e) goes to zero due to (H3|l . 

To get c), one lets C = "^t, 9 = Vt in (12. ip for t G (0,T), and then 
the right-hand member of (12. ip vanishes there. And e) is obtained by 
putting t = in (12. ip and applying a density argument. Finally, d) is 
a consequence of a), e) and c). 
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